Steplike intensity threshold behavior in extreme ionization of laser-driven Xe clusters 
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The generation of highly charged Xe g+ ions up to q — 24 is observed in Xe clusters embedded in 
helium nanodroplets and exposed to intense femtosecond laser pulses (A=800nm). Laser intensity 
resolved measurements show that the high-g ion generation starts at an unexpectedly low threshold 
intensity of about 10 14 W/cm 2 . Above threshold, the Xe ion charge spectrum saturates quickly and 
changes only weakly for higher laser intensities. Good agreement between these observations and a 
molecular dynamics analysis allows us to identify the mechanisms responsible for the highly charged 
ion production and the surprising intensity threshold behavior of the ionization process. 
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Atomic clusters are ideal targets for exploring ultrafast 
laser-driven ionization dynamics of dense matter on the 
nanometer scale Their extreme absorption [2] leads 
to emission of fast electrons d, Q , euv- and x-ray photons 
UlEl, and energetic, highly charged ions While im- 

portant for future particle and radiation sources ^ and 
the generation of plasma waveguides [loj , insight into the 
strong-field ionization mechanisms in laser-cluster inter- 
actions is of great interest also in other fields like laser in- 
duced dielectric modifications and micromachining [Til ] . 

Extreme charging of rare gas clusters in intense 
laser fields has been observed experimentally for over a 
decade [8] . Several effects stemming from the short-lived 
and dense nanoplasma in the Coulomb-exploding clusters 
have been investigated as possible mechanisms for high- 
q ion production. These include tunnel ionization (TI) 
enhanced by fields from neighboring ions, space- charge 



fields, and resonant field amplification [12|, [13|; elec- 



mpliJ 

tron impact ionization (EII) [14[ with atomic ionization 
thresholds lowered by plasma effects in the clusters f]~BI — 
Q3] ; and efficient cluster heating from resonant excitation 
of surface plasmons during cluster expansion [l], [3, with 
possible multiple resonances in core-shell systems [l9[l20|. 
However, the verification of the relative impact of these 
processes, their interrelation and their intensity depen- 
dence is a key outstanding question requiring coordinated 
experimental and theoretical efforts. To date, the aver- 
aging of the ion spectra over the laser focus has impeded 
such a closer connection of experiment and theory. 

We address this long-standing problem by using fo- 
cus z-scanning [2l[ to measure intensity-dependent Xe g+ 
charge spectra resulting from laser excitation of Xe^v 
(N w 150) embedded in helium droplets. The main ex- 
perimental findings are: (i) pronounced intensity thresh- 
olds for strong Xe cluster ionization around 10 14 W/cm 2 , 
(ii) avalanche- like ionization behavior with the appear- 
ance of highly charged Xe g+ with q > 20 even slightly 
above threshold, and (iii) rapid saturation of the charge 
spectra with laser intensity. The intensity window be- 



tween the ionization threshold and saturation of the 
charge spectra decreases with increasing pulse duration, 
leading to a steplike generation of predominantly high-g 
ions (q = 9 — 24, see Fig. [I]) with sufficiently long pulses. 

The ion charge spectra, resolved over a wide range of 
laser intensities, present a unique testing ground for the- 
oretical models. Our molecular dynamics (MD) analysis 
allows us to extract a complete picture of the ionization 
dynamics, backed by good agreement with experiments. 
We reproduce the steplike intensity threshold behavior of 
the Xe charge spectra, and find it is caused by two effects: 
(i) rapid inner ionization of Xe to high-g states through 
an Ell-avalanche sparked by TI of Xe atoms and (ii) sup- 
pression of charge recombination by resonant heating of 
the Xe cluster. The ionization behavior becomes steplike 
for pulse durations being long enough to induce collec- 
tive cluster heating after the avalanche. In addition, the 
MD results show strong ionization of the He shell and, 
due to its rapid expansion, an early plasmon resonance 
(in agreement with Ref. [19|). Though this resonance 
yields high energy capture and further increases Xe inner 
ionization, the Xe cluster dynamics (heating, expansion, 
recombination) is only weakly affected by the He droplet. 

The key to avoiding focal averaging of the ion spectra 
is intensity-selective scanning. We z-scan the laser focus 
with a well-collimated cluster beam, where z is the off- 
set of the intersection point of the optical axis and the 
cluster beam from the focus (Fig. [2k). We thus probe 
a z-dependent intensity distribution [22|, which in prin- 
ciple allows a full retrieval of intensity dependent ion 
yields [23|. For a Gaussian focus profile the peak in- 
tensity on the cluster beam axis is given by the on-axis 
intensity I aK (z) = Iq(1 + z 2 /^ 2 ,) -1 , where Io and 
are the focal peak intensity and the Rayleigh range, re- 
spectively. When reducing the absolute offset \z\, the 
differential interaction volumes dV/dl decrease for all 
present intensities, but increase for the highest one, i.e., 
for / w ^ax(^)- Therefore, z-scanning is very sensitive to 
threshold phenomena making it superior to simply scan- 
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FIG. 1: Xe g+ spectrum from embedded Xejv subject to 1.2 ps 
laser pulses at I ax = 2.5 x 10 15 W/cm 2 showing ions up to 
Xe 24+ and enhanced signal around q = 19. Some charge states 
are compromised by signals from residual gas (as indicated). 



ning laser power, where all interaction volumes are in- 
creased simultaneously. Hence a sharp signal onset at z' 
can be assigned to an intensity threshold at I ax (z f ). 

To produce a stable and collimated molecular beam, 
Xe clusters are grown in He nanodroplets using the pick- 
up method [24[. The droplet beam is produced by su- 
personic expansion of He at 20 bar stagnation pressure 
through a 5 /im nozzle cooled to 9.5 K. In the pick- 
up chamber Xe atoms are successively loaded into the 
droplets. Simulations of the pick-up process [25| indicate 
the formation of XeAr with N up to 500 and TV « 150 
surrounded by ~ 5 x 10 4 He atoms. The cluster beam 
is collimated using a 3 mm orifice, resulting in a similar 
beam diameter in the interaction zone with the laser. A 
30 Hz Ti: Sapphire laser generates 14 m J pulses at 800 nm 
with durations of r = 175 fs - 1.2 ps (FWHM) focused by 
a 40 cm lens (f/33) to 35 /im spot size (zr = 1.2 mm). A 
static 2 kV/cm repeller field accelerates the resulting ions 
into a reflectron time-of- flight mass spectrometer. 

Fig. 1 displays a typical Xe g+ spectrum where charge 
states up to Xe 24+ are resolved. The distribution is 
peaked at high g-values around Xe 19+ and low-q states 
are strongly suppressed. The latter effect is particularly 
pronounced for pulses longer than 1 ps. For z-scanning 
such spectra are measured as function of focus offset, as 
shown for selected pulse durations in Figs.[2j>|2]i. 

Each z-scan profile is a complete image of the r- 
specific, intensity-dependent ion yield up to the focal 
peak intensity. The strong features around z = for 
q < 4 result from TI of Xe atoms in the molecular beam, 
providing an independent measure of the laser parame- 
ters. The asymmetry of the profiles, most pronounced for 
the shortest pulse, is attributed to nonlinear propagation 
effects during laser beam transport through 10 m of air 
which compromise the beam profile and thus the focus 
quality. Irrespective of the asymmetry, the following con- 
clusions can be drawn unambiguously from Fig.[2j (i) In 
all z-scan spectra, Xe 9+ yields from clusters rise sharply 
as function of z with onsets well separated from z = 0. 
For a given pulse length, the signal onset changes only 
weakly with q and roughly coincides with the signal sat- 
uration for Xe + from atomic gas in the beam. This sup- 
ports the claim that TI of Xe atoms in clusters sparks an 
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FIG. 2: (a) Schematics of z-scanning of a Gaussian focus with 
the cluster beam. Contour lines indicate intensity levels rela- 
tive to focus center, (b-d) Measured z-scans of Xe <?+ spectra 
for embedded Xejv (N « 150); pulse durations (r), on-axis 
intensities (lax), and focal peak intensities (Jo) as indicated. 
The strong signals for q < 4 stem from residual Xe gas. 



ionization avalanche producing high-g Xe ions, (ii) The 
maximum q as function of z rapidly saturates, i.e., higher 
on-axis intensities do not yield higher charge states. Most 
notably, the observed threshold intensities for ions up to 
Xe 23+ are near 10 14 W/ cm 2 , whereas the corresponding 
threshold for atomic over-the-barrier ionization exceeds 
10 18 W/cm 2 



Only weak Xe 24+ signal (Figs.Efc, [2i) 
and no ions with q > 24 are found (their generation 
requires 4s electron detachment), which supports that 
atomic shell effects determine the saturation level, (iii) 
With increasing pulse duration the z-scan profiles become 
steplike, i.e., the charge spectra saturate more rapidly. 
This is accompanied by the predominant production of 
high-g ions (cf. Fig. [I]). For the longest pulse, the pro- 
duction of mostly high-g ions sets in promptly and low-q 
ion yields are weak for all z (hat-like feature in Fig.EJi). 

To analyze the ionization behavior and the effect of the 
He shell, we use the MD code from Ref. [I?}- The dynam- 
ics of ions and plasma electrons is described classically 
and linked to an effective quantum mechanical treatment 
of inner ionization for cluster and shell constituents. For 
TI the Ammosov-Delone-Krainov rates [27] employ a lo- 
cal electric field to include plasma field enhancements; 
EII is modeled by an effective Lotz cross section that ac- 
counts for ionization threshold lowering by local fields. 
We use a Gaussian pulse envelope centered at t = 0. 
To determine final ion charge spectra, recombination of 
quasi-free electrons after laser excitation is modeled with 
a simplified scheme: Electrons are treated as recombined 
when localized to an ion after sufficient propagation time. 
More specifically, the inner ionization stage of an ion is 
reduced by the number of electrons that are localized to 
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FIG. 3: Calculated Xe 9+ spectra (color) from pure Xes6i and 
XeseiHesooo subject to 150 fs (top) and 600 fs (bottom) pulses 
vs. laser peak intensity. Curves show averages of final charge 
states (q) (solid) and inner ionization (qa) (dashed). Insets 
depict simulated z-scans (focal peak intensities as indicated). 



it at t rc b = 2 ps, where the ion distribution is frozen. As 
was shown previously, three-body recombination, being 
the dominating process, is efficient only in early stages of 
cluster expansion and can be neglected for t > t vc ^ [17[. 

Fig. [3] shows MD results for Xe56i and XeseiHesooo f° r 
excitation with 150 fs and 600 fs laser pulses. The Xe 
cluster size corresponds to the experimental size distribu- 
tion tail, where highest ionization is expected. For each 
scenario the final Xe g+ distribution is plotted vs. inten- 
sity (color). Curves display average values for inner ion- 
ization ((<fci), dashed) and final charge states after recom- 
bination ((g), solid). Insets contain simulated z-scans re- 
sulting from convolution of the intensity-dependent spec- 
tra with the effective intensity distribution for a Gaussian 
focus using the experimental beam and focus parameters. 

From Fig. [3] we deduce the following: (i) All scenar- 
ios show a sharp onset of cluster ionization that can be 
traced back to TI of Xe atoms, (ii) Beyond threshold, the 
charge distributions shift smoothly to high g-values for 
150 fs, whereas a steplike jump of the spectra to high-g 
ions appears with 600 fs. When comparing the deviations 
of (q) and (qa) it turns out that the impact of recombina- 
tion is responsible for the different saturation behaviors. 
The reduction of ionization stages for 150 fs stretches the 



charge build-up to a finite intensity window, whereas an 
immediate generation of mostly high-g ions and prompt 
saturation of the charge spectra is found with 600 fs. This 
is also clearly visible in the simulated z-scans (insets of 
Figs.[3M3]i). As in the experiment, the z-scan profiles 
show continuous signal from low to high charge states 
for the short pulse, whereas steplike intensity threshold 
features with mostly high-g ions emerge for the longer 
one. (iii) In the presence of the He shell, the saturated 
charge distributions are shifted to higher g-states due to 
more efficient Xe inner ionization. The corresponding z- 
scans (insets Figs. [3b and[3]i) better fit to the experimen- 
tal data (Fig. 2), although the highest predicted charges 
states are about 30% below the experimental values. We 
attribute this to second-order contributions in EII, i.e., 
ionization via intermediate excited states (only the direct 
EII channel is incorporated in the model). 

The reasonable agreement of the MD results with 
the experiment allows us to extract detailed information 
about the various processes involved in high-g ion genera- 
tion. Fig. S] compares the time-evolution of free (dashed) 
and embedded clusters (solid). It reveals the different 
stages of avalanche ionization and the interrelation of TI, 
EII, and heating. The laser intensity (7 x 10 13 W/cm 2 ) 
was chosen such that recombination substantially low- 
ers ionization for the shorter pulse, while the high-g ion 
generation is almost saturated for the longer pulse. 

Irrespective of pulse duration, ionization is sparked 
by TI of Xe atoms. Subsequently, continuous release 
and heating of electrons launch a collisional ionization 
avalanche in the Xe cluster (Figs.Hfci and Hk). When 
present, the He shell is ionized by the avalanche with a 
short delay, in agreement with previous studies [HI, [lo| . 
Interestingly, the evolution of Xe inner ionization is sim- 
ilar with and without helium in early stages. At later 
times, Xe charging is enhanced in the embedded clus- 
ters, which can be traced to resonant collective heating 
of the surrounding He shell. The He resonance occurs 
significantly before the main cluster resonance (the latter 
is present only for r = 600 fs) , as helium expands much 
faster than the heavy Xe ions [19| . The resonant heating 
of the expanded He shell yields ~30% enhancement of 
Xe inner ionization over the free cluster results. 

Most notably, Xe inner ionization ceases as the Xe clus- 
ter expands, i.e., prior to a possible cluster resonance. 
This reflects that high Xe density strongly supports EII 
by higher flux of impinging electrons and stronger ioniza- 
tion threshold lowering [17j. At the cluster resonance, the 
increase of Xe inner ionization is weak (Fig.Hk) despite 
the higher electron temperature in the Xe core (Fig. Eh). 

Whereas the He resonance governs the energy capture 
in the embedded systems, the Xe cluster dynamics is only 
weakly affected by the presence of the He shell. A few 
tens of fs after the He resonance, the electron tempera- 
ture in the Xe core evolves very close to the free cluster 
results (Figs.|4]i and Eh). This shows that the long-term 
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observed. Our MD results indicate that avalanche clus- 
ter inner ionization proceeds before a cluster plasmon 
resonance may occur and is enhanced by collective heat- 
ing of the surrounding He shell. The main role of the 
Xe cluster resonance is suppression of charge recombina- 
tion, i.e., conservation of the ion distribution produced by 
avalanching. With higher focus quality in future experi- 
ments, z-scanning promises access to the full intensity- 
resolved ionization rates of clusters or other complex 
nanosy stems. The use of shorter pulses would offer in- 
sights into the transition between avalanche and multi- 
photon driven material damage, a topic of great interest 
in micromachining and laser processing of dielectrics. 

Main parts of the helium droplet machine have been 
provided by J. P. Toennies (MPI Gottingen). Support by 
the DFG within the SFB 652 and by the HLRN Com- 
puting Center are gratefully acknowledged. 



FIG. 4: Evolution of inner ionization, absorption, clus- 
ter/matrix radii, and average Xe-core electron energy (top 
to bottom) for free and embedded Xes6i in 150 fs and 600 fs 
pulses at Io=7 x 10 13 W/cm 2 . Panels (a) and (e) also show 
the effect of recombination for Xe, evaluated at t rc b = 2 ps. 



electron temperature evolution in the Xe core is mainly 
determined by direct laser heating of the Xe cluster, i.e., 
mostly unaffected by the He shell. As the rate for three- 
body recombination is strongly reduced at high temper- 
ature [111, resonant heating of the Xe cluster leads to 
conservation of high-g ions (Fig.|4^). In contrast, the ion 
charge states are substantially reduced by recombination 
in the absence of a cluster resonance (Fig.|4^i). Hence, 
resonant excitation of the cluster plasmon is the key to 
avoiding a strong reduction of the high charge states by 
recombination. The main effect of the He shell is an inner 
ionization enhancement in early stages of the interaction. 

The interplay of avalanching and recombination deter- 
mines the threshold behavior. With short pulses, Xe 9+ 
spectra rise smoothly with intensity above threshold un- 
til avalanching begins so early in the pulse that reso- 
nant heating in the trailing edge becomes efficient. With 
sufficiently long pulses (r > 1 ps), suppressed recombina- 
tion due to resonant heating occurs whenever avalanching 
takes place (for I>10 14 W/cm 2 ), i.e., the threshold be- 
havior becomes steplike. Higher intensities do not change 
the ion spectra significantly, as avalanching and resonant 
heating just proceed earlier in the pulse. Only at much 
higher intensity, atomic TI is expected to further increase 
charge states, but requiring / > 10 18 W/cm 2 for q > 20. 

In conclusion, we used z-scanning to measure intensity- 
resolved ion spectra from embedded Xe clusters. Inten- 
sity thresholds for producing Xe 9+ ions with q > 20 of 
only 10 14 W/cm 2 and a steplike ionization behavior are 
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